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Abstract Hexagonal BaTiC>3 undergoes a structural phase transition to 
an orthorhombic CTX1\ phase at To = 222 K. The transition is driven 
by a soft optical mode with Ei u symmetry whose couplings force the ap- 
pearance of a spontaneous E2 g strain (improper ferroelastic character). 
Staying within the same Em subspace, the system could in principle set- 
tle into a second (not observed) orthorhombic phase (Cmc2\). We have 
carried out a first-principles investigation of these questions, studying 
the structure of the existing C222i and the virtual Cmc2\ phases, and 
describing the spontaneous E<i g strain in accord with the experimental 
observations. In addition, we show that the occurrence of C222i instead 
of Cmc2\ cannot be explained by the E^ u soft modes themselves and, 
therefore, must be related to their couplings with secondary order param- 
eters. A more detailed analysis proves that the Em strains do not account 
for the experimental preference. 
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INTRODUCTION 



Barium titanate has two structural polymorphs: the well known per- 
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ovskite type (p-BT) and its hexagonal modification (h-BT) with 30 
atoms per unit cell. h-BT Undergoes two zone-center structural 
phase transitions: at To = 222 K from the hexagonal P&^/mmc 
phase (I) to an orthorhombic C222i (II), and at Tc = 74 K to a 
monoclinic P2\ phase (III) . @ The transition at To presents two spe- 
cially interesting features: 

• It is driven by a "silent" soft optical mode of E 2u symmetry, 6 cou- 
pled quadratically to the E 2g acoustic modes (strains) rji — r\ 2 and 
t]q, which spontaneously acquire non-zero values at the transition.® 
In the usual terminology, the transition is then improper ferroelas- 
tic. The coupling of the E 2g modes with the soft optical modes has 
been studied spectroscopically.'Q'ti 

• Depending on the direction taken by the bidimensional E 2u order 
parameter, the system can find itself in structures of three different 
symmetries: the experimentally observed C222! and P2 1; as well 
as a second orthorhombic Cmdl\ phase (II') . The preference of II 
over II' must ultimately be dictated by the energetics of the E 2u 
instability and by its couplings with secondary order parameters. 

In this paper we present a summary of the results of our investigation 
into the E 2u optical instability in h-BT, including its couplings with 
the spontaneous E 2g strains, and into the issue of the relative stabil- 
ity of the two accessible orthorhombic phases, using first-principles 
density-functional methods. To our knowledge, this is the first ab 
initio study of hexagonal BaTiCV 



TECHNICAL DETAILS 



For our calculations we have used density-functional theory within 
the pseudopotential approach and the local-density approximation. ® 
We have employed Vanderbilt's ultrasoft pseudopotentials,® which 
are ideally suited to deal with first-row elements and transition met- 
als. The electron wave functions were expanded in a plane wave 
basis with an energy cut-off of 25 Rydberg, and the Brillouin zone 
sums were calculated by a 3 x 3 x (2(+0.5)) Monkhorst-Pack spe- 
cial k-point mesh.il We checked the validity of these approximations 
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by a convergence analysis. This first-principles framework has been 
successfully used for the investigation of structural properties in per- 
ovskite oxides, o 



RESULTS AND DISCUSSION 



In order to proceed, the reference hexagonal structure must be deter- 
mined. In I the occupied Wyckoff positions are fixed by symmetry up 



to five free parameters. We used the systematic procedure of Ref. |10 
to relax these five A lg optical modes so that the forces on atoms were 
below 10~ 4 eV-bohr^ 1 . Our results for the fractional atomic coor- 
dinates agree with experiment0 to within 0.6%. Due to the high 
computational cost of a complete analysis, we used the experimental 
values for the unit cell parameters (a = 10.77 bohr and c/a = 2.456). 

In what follows we work with a symmetry-adapted Taylor expansion 
of the potential energy of the system around the reference structure, 
as a function of the above mentioned relevant variables: the E 2u soft 
optical and the E 2g acoustic modes. It is convenient to decompose 
this expansion into: 

E E\ -^optical -^strain -^coupling; (1) 

where E\ is the energy of the reference structure I. 

As a pre-requisite for the study of -E op ticai, we needed to check that 
the ab initio methods predict the experimentally observed instability. 
We solved the associated eigenvalue problem (there are seven zone 
center E 2u symmetry modes) and indeed found only one instability 
with spring constant —0.743 eV-bohr~ 2 . (We consider the Fourier 
transform of the so called force-constant-matrix and not the dynam- 
ical matrix.) The corresponding (degenerate) eigenvectors determine 
the unstable two-dimensional subspace. Using as coordinates the re- 
spective amplitudes in polar form (Q 1 =: rcosp and Q 2 ='■ r sirup), 
we can write: 

Optical = \A 2 r 2 + + \ [As + A' Q cos{Qp)}r & , (2) 
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where <p = {0, ±7r/3, ±2n/3, tt} and <p = {±7r/6, ±7r/2, ±57r/6} cor- 
respond to II and II' respectively, while for a general ip the system 
is in III. Eq.(||) shows that the sign of A' 6 determines the relative 
stability of the orthorhombic minima (with A 2 < 0)|] 
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FIGURE 1: Potential energy well of the E 2u instability of I, in the 
(p = direction (corresponding to II). The ab initio data is fitted to 
a quartic and Ej is set to zero. 

Figure ([!]) shows -E op ticai as a function of r with ip = 0, the case of 
II. (The data for II' are identical at the scale of the figure). In order 
to model these data, two facts must be taken into account. First, 
quadratic (v4 2 ) and quartic (A4) terms are symmetry-forced to be 
the same in both phases. Second, the fitted A 2 must be in good 
accord with the previously obtained eigenvalue. We determined: 
A 2 = -0.72 ± 0.02 eV-bohr" 2 and A 4 = 2.56 ± 0.07 eV-bohr" 4 , with 
the energy minima in |r m j n | = 0.52 ± 0.01 bohr [E m in = — 50 ± 4 
meV). We could obtain no convincing fitting for the 0(r 6 ) term and 
thus concluded that the quartic approximation constitutes the best 
model to describe the E 2u instabilities, i.e., the soft modes themselves 
cannot account for the differences between the II and II' minima, so 
their couplings with secondary order parameters must be responsible 
for the experimental observation of II. 

This calculation constitutes a prediction of the structural parame- 
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■f For the energetic analysis of the monoclinic phases, higher order terms 
should be considered. 
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ters of the orthorhombic Clll^ phase, which to our knowledge has 
not yet been determined experimentally. (Further details about the 
structure will be published elsewhere.) 

Turning to the pure elastic energy term in Eq. (12) , it is convenient to 
use polar coordinates also for the E 2g strains (771 — r/ 2 =: pcoscj) and 
r/ 6 =: psintfi), and write: 

Strain = <V + C 3 COs(3(f>) P \ (3) 

For a general the system is distorted to III (due to rje), while for 
= {0, 7r} it goes to Cmcm, which is a supergroup of both II and 
II', so 771 — r]2 is the only secondary order parameter we take into 
account. We studied the response of the system to this strain with 
the optical phonon amplitudes set to zero (and imposing 771 +r? 2 = 0, 
so that only deformations with E 2g symmetry are considered), and 
concluded that the best model reduces to C 2 = 130 ± 2 eV. (The 
0(p 3 ) term cannot be reliably fitted.) 

Finally, we examined the coupling between the E 2u soft optical and 
the E 2g acoustic degrees of freedom, considering: 



p - pW _L p( 2 ) 

-^coupling -^coupling ' "^coupling 

Kipling = a 12 cos{2ip + 0)pr 2 

Kipling = {aiACOs{2if + 0) + (3 u cos{Aif + 0)} pr 4 

+ {a 22 + (3 22 cos [2(<p - <f))}}p 2 r 2 

+ {"24 + P24COS [2(<p - 0)] + 724COS [2(2<p - 0)]} p 2 r 4 , 

(4) 

The first term -^coupling * s the simplest form of coupling usually con- 



sidered in the literature, O'EJ and indeed it provides us with a very 
good fit of the ab-initio energies in the presence of E 2g distortions 
(for both the II and II' phases). Figure @ shows the main re- 
sult of including this minimal strain coupling into the energy anal- 
ysis: the minima shift their positions and get deeper. As table 
(TJ) shows, the configurations of minimum energy are located at 
\ r min\ = 0.55 ± 0.01 bohr, and the associated spontaneous strain 
is Pmin = (5.6 ± 0.2) x 10~ 3 with = and n for II and II' respec- 
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tively. We found no experimental measurement of the aci 2 parameter; 
nevertheless, its sign is correct as compared with results in Ref. [TJ. 



The location of the wells provides us with a more complete (this 
time including strain) prediction of the structure of the orthorhom- 
bic phase. 
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FIGURE 2: (r, p) Energy map of II. We show the equi-energy lines 
of a model in which only the A 2 , A 4 and C 2 terms (with the values 
of the text) are present (dashed) and those corresponding to the fit 
in table (|I|) (solid). Energies are given in meV and E\ is set to zero. 

However, as the well depths are within this minimal model identi- 
cal for both II and II', we still have found no reason for Nature's 
preference for the first. To further analyze this issue, we included 
-^coupling i n Eq.(|J) in our fitting procedure and found that only the 
0(p 2 r 4 ) terms improve the model. We got a 2 ± + /?24 + 724 = 7.3 ± 0.2 
eV-bohr" 4 for II and a 2i - (3 2 i + 724 = 7.9 ± 0.2 eV-bohr~ 4 for II', 
but the difference is so small that it produces no definite distinction 
between the two minima. Thus, it can be concluded that neither 
the E 2u instability itself nor its coupling with E 2g strains justify the 
existence of II instead of II'. 



CONCLUSIONS 



We have presented a first-principles study on the stability of hexag- 
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ii ir ii & ir 


C 2 (eV) 
«i2 (eV-bohr -2 ) 


134 ±2 133 ± 2 133 ± 2 
-5.0 ±0.1 -4.8 ±0.1 -4.9 ±0.1 


pmin (10" 3 ) 
\r min \ (bohr) 
E min (meV) 


5.6 ±0.2 -5.4 ±0.2 (±)5.6 ± 0.2 
0.55 ±0.01 0.55 ±0.01 0.55 ± 0.01 
-54 ±4 -54 ± 4 -54 ± 4 



TABLE 1: Best fits for Eq.(|T]) when only a minimal coupling (a\ 2 
term) is considered in the model. We present the fits corresponding 
to II and II' separately (so that the fulfillment of symmetry restric- 
tions can be observed), plus a final result (II & II') including aggre- 
gate data for the two phases. In all cases, parameters A 2 and 
were fixed to the values in the text. The quality of the fits is highly 
sensitive to small variations in C 2 , so we left it free (the obtained 
results are within the error of our initial estimation). 

onal barium titanate when distorted by the zone center E 2u soft 
optical and the E 2g acoustic modes. 

We have shown that first-principles calculations are able to predict 
the experimentally observed zone center E 2u instability and to define 
the corresponding unstable subspace of distortions. We have exam- 
ined the anharmonicity of the soft modes and obtained a prediction of 
the structure of the observed C222i and of the virtual Cmc2\ phases, 
including the appearance of a spontaneous E 2g deformation which is 
coupled to the soft mode. Moreover, we find that the two accessi- 
ble orthorhombic minima are not definitely distinguished by the E 2u 
instability itself, which points to the importance of secondary order 
parameters to decide this issue. However, the effect of the coupled 
E 2g strains is not large enough to explain the experimental obser- 
vation of C222i instead of Cmc2i, which must then be due to the 
couplings of the soft mode with other secondary order parameters. 
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